Electron transport in disordered graphene 
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We study electron transport properties of a monoatomic graphite layer (graphene) with different 
types of disorder. We show that the transport properties of the system depend strongly on the 
character of disorder. Away from the half filling, the concentration dependence of conductivity is 
linear in the case of strong scatterers, in line with recent experimental observations, and logarithmic 
for weak scatterers. At half filling the conductivity is of the order of e 2 /h if the randomness preserves 
one of the chiral symmetries of the clean Hamiltonian; otherwise, the conductivity is strongly affected 
by localization effects. 
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I. INTRODUCTION 

Recently, Novoselov et al have succeeded in fabrica- 
tion of monolayer graphite (graphene) samples.- 1 Subse- 
quent transport measurements^*^^ have shown that 
graphene is a conductor with remarkable electronic prop- 
erties. These experimental discoveries have triggered 
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crs in graphene have a relativistic (Dirac) spectrum^2ii£ 
which makes the transport properties of this material 
highly interesting from the point of view of both fun- 
damental physics and potential applications. It is widely 
believed that graphene-based devices may be of outstand- 
ing importance for future nanoelectronics. 

This work has been motivated by the following two ex- 
perimental observations!^ First, it was found that the 
graphene conductivity is linear in the concentration of 
carriers (counted from the half filling) with a high accu- 
racy. Second, it was found that at half filling the conduc- 
tivity (per spin direction and per valley) is close to e 2 /h 
and does not show any definite temperature dependence 
in a broad temperature range. The aim of this paper is 
to analyze what one should expect for conductivity from 
the theoretical point of view and whether these theo- 
retical predictions may be compatible with experimental 
findings. We will see that, in view of the unconventional 
character of the graphene spectrum, the theoretical re- 
sults depend crucially on the nature of disorder. 

The structure of the paper is as follows. In Sec. ITT1 we 
introduce the model describing electronic properties of 
graphene with various types of disorder. In Sec. IIIII we 
analyze the dependence of conductivity on the electron 
concentration away from the half-filling point. We con- 
sider the two limits of weak and strong scatterers and 
construct the corresponding "phase diagram" . Section 
FVlis devoted to the conductivity at half filling under the 
assumption that the disorder preserves one of the chiral 
symmetries of the Dirac Hamiltonian. Our findings are 
summarized in Sec. [V] Some technical details are pre- 



sented in two Appendices. 



II. THE MODEL 
A. Clean graphene 

The carbon atoms of graphene are arranged in the hon- 
eycomb lattice (see Fig. QJi) with the period a = 2.46 A. 
Each carbon atom of intrinsic graphene has one valence 
electron forming the 7r-bonds to the three neighbors. The 
electronic spectrum of graphene is well described by the 
tight-binding modelSS taking into account the nearest- 
neighbor hopping. The first Brillouin zone for this system 
has a form of a hexagon (see Fig. QJ)) with the distance 
fco = 2h/3a from the center to the apex. The honeycomb 
lattice contains two sites per elementary cell. This per- 
mits the grouping of all the atoms into two sublattices, 
A and B. The nearest neighbors of an atom from the 
sublattice A belong to the sublattice B and vice versa. 
The symmetry group of the honeycomb lattice contains 
an element swapping the two sublattices. Hence, for each 
value of quasimomentum k within the Brillouin zone, two 
states exist with the energies ±£'(k). These two spec- 
trum branches are degenerate at the isolated points in 
the corners of the Brillouin zone, E(ko) = 0. With one 
electron per site the system is exactly in the half-filling 
state when the nodal points of the spectrum lie at the 
Fermi level. Among six apices of the hexagonal Brillouin 
zone only two are nonequivalent. They are referred to 
as K and K 1 . The electrons with momentum close to 
these two points, and hence with low energy, are rele- 
vant in studying the physics of the system for electron 
concentrations not too far from the half filling. 

The tight-binding Hamiltonian is a 4 x 4 matrix oper- 
ating in the AB space of the two sublattices and in the 
K-K 1 space of the valleys. Therefore we introduce the 
four-component wave function 



* = {4>ak, <I>bk, 4>bk>Aak'Y 



(1) 
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FIG. 1: (a) Honeycomb lattice of carbon atoms of graphene. 
Solid and open circles denote the atoms of A and B sublattices 
respectively, (b) The first Brillouin zone of graphene. The 
nodal points of the spectrum are located in the corners of the 
zone. The two nonequivalent nodal points are denoted as K 
and K'. 



In this representation the Hamiltonian has the form 

H = V Q T 3 (Tk. (2) 

Here t 3 is the third Pauli matrix in the K-K 1 space and 
er = {(Ti, (T2 } is the two-dimensional vector of Pauli ma- 
trices in the AB space. The Fermi velocity in graphene is 
vo ~ 10 8 cm/s. In fact, the form of the Hamiltonian 
is universal and does not rely on the tight-binding ap- 
proximation. The degeneracy of the spectrum in K and 
K' points is provided by the two-dimensional representa- 
tion of the honeycomb lattice symmetry group while the 
expression J5J is the first-order fc-expansion near these 
points. As k is increased the higher-order nonuniversal 
terms of this expansion come into play. For our pur- 
poses, it will be sufficient to introduce the high energy 
cut-off A and to assume the spectrum to be linear up to 
|k| = A/vq. Indeed, all divergent momentum integrals 
appearing below have the logarithmic character; thus, 
details of the high-energy regularization are irrelevant. 
The Green function for the Hamiltonian J5J of the clean 
graphene reads 



E + t>oT 3 crk: 
(e±i0) 2 -v 2 k 2 ' 



(3) 



B. Potential disorder 

We incorporate now disorder in the model. Let us 
consider first the impurities modifying the potential on 
nearby lattice sites. A detailed description due to Mc- 
Cann and Fal'ko^ contains 10 real parameters for the po- 
tential of a single impurity. In the present paper we will 
use the simplified model introduced by Shon and Ando in 
Ref. 0, which retains the essential physics of the prob- 
lem. This model treats impurities in the framework of 
the same tight-binding approximation as was used for 



the pure system. An impurity is placed at a site of the 
lattice and has a potential U (r) . We use the two discrete 
Fourier transforms of this function with respect to the 
two sublattices. 
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-zqr 



(4) 
(5) 



The summation runs over all elementary cells of the hon- 
eycomb lattice, and the vector m points from the A sub- 
lattice site to the B sublattice site of the same elementary 
cell (see Fig. QJi). The quantity Uq is the scattering am- 
plitude for the electrons of the same sublattice where the 
impurity resides, while ?7 q is the scattering amplitude for 
the electrons of the other sublattice. 

Assuming C/ q and C/ q are slow functions of the mo- 
mentum q, we keep only two values of these amplitudes 
for intravalley, Uq and Uq, and intervalley, U ko and U' k , 
scattering. The hexagonal symmetry of the honeycomb 
lattice makes the amplitude U' k vanish while the three 
other amplitudes are real. Thus we are left with the three 
parameters of an impurity potential. This is straight- 
forward to put them in a matrix in the 4-dimensional 
representation Q. If the impurity site belongs to the 
sublattice A and to the elementary cell rj, the scattering 
matrix takes the form 
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For the impurity located in the sublattice B, we have 



(6) 
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U U ko e~ 2lk ^ 
U ko e 2 ^ Uo 


(7) 

If the potential disorder is weak and obeys Gaus- 
sian distribution, the only relevant quantity is the auto- 
correlation function of the second order (Vq ® V_ q ). We 
denote the impurity concentration by ni mp and obtain af- 
ter averaging with respect to positions of the impurities 



(Vq «) V- 



'imp /i/A 



(F q A (r 4 ) ® y A q (rO + V*(ri) ® Vl B q (rO) 

= 27TU 2 ja <7 Q Tq (g) Cr T + 7 2 <7 3 T 3 <g> (T3T3 

#1 r 

+ — L<Tm ® (Tin + 01T2 ® C1T2 

+ CT2T1 (g) fT 2 Tl + C2T2 <g> (T2T2] |. (8) 

Here we introduce the following three dimensionless pa- 
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rameters 



a 



"'imp 
Tl'imp 
^imp 



(Uo Utf 



(9a) 
(9b) 
(9c) 



While the notations in Eq. may seem strange at 
this stage, they will be explained later when we consider 
randomness of a broader class. To further simplify the 
calculations, we will concentrate on the two limiting cases 
of short- and long-range potential disorder4i 

The short-range impurity scatters electrons in the 
same sublattice only. It is equivalent to a potential 
shift at a particular lattice site. The amplitudes are 
U = U ka = U/2 and Uq = 0. Thus we are left with 
a single parameter U. The parameters Eqs. © obey the 
relation ao — j z — (3±/2. 

The long-range impurity scatters electrons in both sub- 
lattices equally but only within one valley. The scatter- 
ing length is large in comparison with the lattice constant 
but is still smaller than the Fermi wavelength. The am- 
plitudes are Uq = Uq = U, Uk = 0. We have again the 
single parameter U as in the case of short-range disor- 
der. Among the parameters @ only «o is not zero in 
this case. 



C. Generic disorder and chiral symmetries 

Let us turn now to the analysis of the symmetries of 
the clean graphene Hamiltonian @. First, the system is 
obviously uniform and isotropic. Any disorder considered 
in this paper preserves these symmetries on average, so 
we do not pay much attention to them here. Second, due 
to the two valley structure of the electron spectrum the 
whole SU(2) symmetry group exists in an isospin space 
of the valleys. The generators of these group are2& 



A, 



03 Tl, 



(10) 



These three operators commute with the Hamiltonian 
and anticommute with each other. There are other three 



matrices ^ x ,y,z introduced in Ref. |26 



02T3, 



£z = a 3 T . 



(11) 



These operators generate an additional SU(2) group of a 
pseudospin. They do not commute with the Hamiltonian 
PJl , however any of these matrices commute with any of 

Third, the time inversion operation (we denote it To) 
in the representation Eq. reads 



A I ► G\T\A 0~\T\. 



(12a) 



The Hamiltonian Eq. J2J is invariant under time inversion 
(note that the momentum operator changes sign under 



transposition). Combining the To operation with any of 
A x ,j/,2 from Eq. I|1U|) we produce three additional symme- 
try operations 



A h-> o- 2 t A 7 '02T0, 
A h-> o- 2 t 3 A t o 2 t 3 , 

A I— > G\T 2 A T G\T 2 . 



(12b) 
(12c) 
(12d) 
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Finally, there is one more — namely, chiral — symme- 
try Co, and its three counterparts generated by simulta- 
neous application of Co and A- x<VtZ : 



(13a) 
(13b) 
(13c) 
(13d) 



The chiral symmetry Co can be viewed as the basic 
chiral symmetry of the Hamiltonian @. Indeed, Co is 
distinguished by the fact that it is directly produced by 
the Hamiltonian © as ia 3 TQ — VQ 2 (dH/dk x )(dH/dk y ), 
while other chiral symmetries require a rotation in the 
isospin space. 

Generally the chiral symmetry implies that the Hamil- 
tonian takes block-off-diagonal form under a proper uni- 
tary transformation. A generic disorder preserving C z 
symmetry can have only off-diagonal matrix elements in 
the AB space of sublattices. Some specific examples of 
chiral symmetry are (i) bond disorder due to distortions 
of the lattice (C z symmetry) , (ii) random magnetic field 
(all four symmetries Co, x ,y jZ ), (hi) dislocations, that are 
equivalent to a random non-Abelian gauge fielo»^2i (Co), 
(iv) infinitely strong short-range on-site impurities (C z ). 
In the latter case an electron cannot occupy the impurity 
site, implying that all the bonds adjacent to the impu- 
rity are effectively cut. Any potential disorder other than 
the described extreme case violates all chiral symmetries. 
The symmetry is also broken by a non-zero chemical po- 
tential. Thus, the impact of the chiral character of dis- 
order will be particularly important at the degeneracy 
point e = 0. In Sec. II VI we consider various effects of chi- 
ral symmetry on the density of states and conductivity 
of graphene. 

The average isotropy of the disordered graphene im- 
plies that A x and A y symmetries of the Hamiltonian 
are present or absent simultaneously. Below we com- 
bine them into a single notation Ax, and proceed in the 
same way with T± and C±. In Table [I] we list all possible 
matrix structures of the disorder [in the representation 
defined by Eq. Q] along with their symmetries. There 
are 9 different structures altogether iSi Those 5 of them 
that do not violate time inversion symmetry coincide^ 
with ones considered by Aleiner and Efetov^S We also 
give the notations of Ref. 4{J where the disordered Dirac 
Hamiltonian obeying C z chiral symmetry was considered. 
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TABLE I: The symmetries of various disorders in graphene. The first five rows of the table contain disorders preserving time 
inversion symmetry. They were considered in Ref. 1.3(1 Next four rows are occupied by disorders violating time inversion 
symmetry. We present the matrix structure of the disorder in two forms: by matrices <JiTj and by matrices EiAj as in Ref. 12^. 
The notations we use for the amplitudes of the disorder in Gaussian limit are listed in the third column. The letters a, f3, and 7 
correspond to Ao, A x ,y, and A z components of the disorder Hamiltonian respectively, while the subscripts 0, _L, and z indicate 
the structure in £0, E I|S , and E z domain. In the fourth and fifth column we give alternative notations from Refs. [3(] andEljl 
Our notations are close to those of Ref. l3*cl the only difference is in the case of a fully diagonal potential: our parameter ao 
corresponds to 70 from Ref. I30L while we use 70 for the disorder oqtz discriminating the two valleys. 
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III. CONDUCTIVITY FAR FROM THE 
DEGENERACY POINT 



In this section, we will study the concentration depen- 
dence of the conductivity far from half filling, when the 
size of Fermi circles around K and K' points is large in 
comparison with the inverse mean free path. The dimen- 
sionless Drude conductivity (measured in units of e 2 /h) 
is then large^i so that, at realistic temperatures, one 
can neglect as a first approximation the quantum cor- 
rections related to localization. As a starting point, we 
will employ the self- consistent T-matrix approximation 
(SCTMA 45 ) that takes into account all orders of scatter- 
ing at an impurity. It will allow us to study the whole 
"phase diagram" including the limits of weak (Born) and 
strong (unitary) scatterers and the crossover between 
them. We will discuss the status of SCTMA in Sec. llTTBl 
where we will show that while it is not quantitatively 
justified in the Born regime, it yields a qualitatively cor- 
rect behavior of the conductivity far from the degeneracy 
point. 



A. Self-consistent T-matrix approximation 

1. Potential disorder 

Let us first consider the disorder induced by randomly 
located impurities which create the potential given by 
Eqs. 10 and J7J. The sum of all scattering orders deter- 
mines the complete impurity's T matrix, as represented 
graphically in Fig. [3 Averaging the T matrix with re- 



spect to the position of the impurity, we find 
2% 



TO) 



1 - %9 
U + U ko 



U - U ko 



l-(U + U ko )g 1-(U -U ko )g 
with g being the integral of the Green function, 

d 2 k 



9(e) 



:G(e,k). 



(14) 



(15) 



This quantity has trivial matrix structure due to the 
angular integration. The electron's self energy is deter- 
mined by the average value of the T-matrix Eq. (|14H 



where 



S(£)=n imp (T( e )) ) 
is the concentration of impurities. 



Inserting (|16fl into the bare Green function 
£ — S(e) + voT 3 crk 



G(e,k) 



my 



v^k 2 



(16) 



(17) 



and calculating the momentum integral in Eq. i|15fl . we 

get 



9(e) 



e-S(e) , 



-A 2 



E(e) 



(18) 



The logarithmic divergence is cut at the momentum 
A/u . The sign of the imaginary part of <?(e), and hence 
of the self energy, is determined by the type of the Green 
function (advanced or retarded) we are considering. 



+ 



FIG. 2: Graphical representation of the T matrix describing 
the electron scattering off an impurity. 



The equations 1)14)) . (|16[) . and l|18l) form a closed set 
that self-consistently determines the self energy E(e). 
These equations take into account all the diagrams with 
non-intersecting impurity lines. In the two extreme cases 
of short- and long-range potential impurities the self- 
consistency equation reduces to the form 



S( £ ) 



A[l-Ug(e)] 



long-range; 



short-range. 



(19) 



Once these equations are solved, one can find the den- 
sity of states and the conductivity of the system. The 
density of states (per one spin component) is 



1 



ImTr 



d 2 k 



G R (e,k) 



Img R (e). 



(20) 

The conductivity at zero frequency and wave vector is 
given by the Kubo formula, 



a a ' j (e) = - / cf (r-r') 



Tr( j a Im G R (e; r, r')/ Im G R (e; r', r) }. (21 



Due to the linear dependence of the Hamiltonian on k, 
the current operator is independent of the momentum: 



dH 

J = e— = ev r 3 a. 



(22) 



This results in the absence of the diamagnetic term in 
the expression for conductivity. 

Equation l(2*T)l includes averages of the type (jG R jG A ) 
along with (jG R jG R ) and (jG A jG A ). The first one is 
large in the metallic regime when the energy e is far 
from the degeneracy point, while the two others give a 
contribution of the order of conductance quantum e /ft. 
Therefore we will neglect those two in this section. 

As discussed above, we will use the Drude approxi- 
mation for the conductivity, neglecting weak localization 
corrections. Graphically this is equivalent to the sum- 
mation of the diagrams shown in Fig. [3J Due to the 
vector nature of the current operator in the vertex, only 
the diagonal parts of the Green functions contribute to 
the result. We introduce the notation 



U RA (e) 



d 2 k 
(2n) 2 



diagG fl (e,k)diagG" 4 (£,k). (23) 



FIG. 3: Diagrams for the Drude conductivity including the 
vertex correction. 



The sum of the ladder diagrams in Fig. |3 give the correc- 
tion to the current vertex. We will use a special notation 
V for this vertex correction factor. In the limit of the 
short-range potential disorder, we have V = 1. In the 
opposite long-range case, the summation of ladder dia- 
grams yields 

' (24) 



V 



1 - 



^ A 



The resulting Drude conductivity has the form 

4 



o-(e) = -eVVn 

7T 



ha 



(25) 



In the following sections we will solve the self-consistency 
equations and find the density of states and the conduc- 
tivity in various limits. We will also analyze the correc- 
tions to the SCTMA coming from diagrams with inter- 
secting disorder lines. 



2. Generic Gaussian disorder 

In the case of generic weak disorder, one can use a more 
general equation for the self energy taking into account 
all possible disorder amplitudes (listed in Table UJ) in the 
framework of the Born approximation, 



E(e) = 2nv 2 a g(e). 



(26) 



Here a is the total strength of the disorder, i.e. the sum 
of all amplitudes from Tabled 

a = a + (3 + 70 + a± + Px + 7-L + a z + Pz + 7z- (27) 

This quantity is relevant for thermodynamic properties of 
the system with Gaussian disorder. The vertex correction 
is given by 

V= l-47rug(a-atr)n^' (28) 

where 

1 3 

atr = ^ ( a o + A) + 7o) + a± + (3± + 7_l + - (a z + f3 z + j z ) . 

(29) 

As will be seen below [Eq. (|37|l ] , atr governs the transport 
properties of the Gaussian disordered system. 

Using Eqs. ©, 123, and we get for weak 

long-range and short-range potential disorder considered 
above in Sec. IIII A II 



a = 2a t 



a = atr 



n- U 2 

_ '"amp ^ 

= 2irv 2 ' 

^imp U 

8ttv 2 ' 



long-range; (30) 
short-range. (31) 
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B. Born limit: Weak scatterers 



1. Self-consistent Born approximation 





(c) 



The simplest situation is the Born limit of weak scat- 
tering. Only the lowest scattering order is relevant in 
this case. We expand the right-hand side of the equation 
()19JI to the second order in the scattering amplitude U. 
(The first-order term is a constant which is absorbed in 
renormalization of the chemical potential.) This yields 
the self- consistent Born approximation (SCBA). 41,54,55 
In this limit, we deal with a generic weak disorder de- 
scribed by all 9 parameters listed in Table [I] The results 
for the special case of potential disorder are easily re- 
stored from the general results with the help of Eqs. 

The SCBA equation has the form 



E( £ ) 



a r 



E(e)] log 



-A 2 



(32) 



This e quat ion was studied numerically by Shon and Ando 
in Ref. l41t we treat it below by analytical means. 

Weak disorder introduces an exponentially small en- 
ergy scale 



To = Ae 



-l/c 



(33) 



At large energies, e 3> T , we solve the equation 
by iterations, while at low energies the solution is found 
in the form of a series in powers of e. The resulting self 
energy is 5 ^ (upper sign - retarded, lower sign - advanced) 



FIG. 4: Diagrams yielding logarithmic corrections to the con- 
ductivity not included in the SCBA. 



With the help of Eqs. lt2"5|l and (|2"4"jl we find the general 
expression for the SCBA conductivity 



cscba(£) = — o 

7T Z 



atr £ — a Re E(e) 



(37) 



Here the first term comes from the retarded-advanced 
(RA) sector whereas the second term (unity) is the con- 
tribution of RR and AA correlators. At high energies 
(e ^> To) the SCBA conductivity is governed by the RA- 
term and takes the form 



o~scba{s) 



it 2 



I 



a 2 , A 

— 1o Sm 
atr £ 



(38) 



The found conductivity shows a logarithmic energy de- 
pendence above an exponentially small energy scale. At 
the half filling, e <C r , the SCBA yields oscba = 
2e 2 /n 2 h. This value of conductivity includes contribu- 
tion of the form (jG R jG R ) and (jG A jG A ), which were 
discarded at e ^> Tq . A conductivity value of the order 
of e 2 /h does not make much sense in the present context, 
in view of the localization effects. We will return to this 
issue in Sec. IIVI 



E(e) 



, A i7ra|e| 

" a£l0g R T ^ 



\e\ <r ; 



1 + 2a log 



A 



(34) 



Substituting Eq. Ip51)l into Eq. lj2T))l. we get the density of 
states, 



Pscba(e) 



2|lmE(e)| 
ir 2 v 2 a 

2Tn 



= < 



2 2 5 



A 

2a log — 

£ 



£< r ; 



e> r . 



(35) 



At high energies the found density of states is close to its 
value in clean graphene, po(e) — \e\/irv 2 . 

To evaluate the SCBA conductivity, we first find the 
polarization operator (|23|l . 



Il RA (e) 



1 



Anv 2 a e — Re E(e) ' 



(36) 



2. Logarithmic corrections and renormalization group 

The leading term in the Drude conductivity (|37|l is 
proportional to a^ 1 and is independent of energy. The 
SCBA gives also logarithmic corrections, which are small 
at large energies. There exist, however, other contribu- 
tions of the same order that are not included in the SCBA 
calculation, see Fig. 21 An efficient tool for resummation 
of the logarithmic contributions in all orders is the renor- 
malization group (RG). For the case of 2D Dirac fermions 
subjected to various types of disorder it was developed 
by Dotsenko and Dotsenko 5 ^ for the random bond Ising 
model, by Ludwig et a& in the context of the quantum 
Hall effect, by Nersesyan et aJSi and Bocquet et a& in 
application to dirty superconductors with unconventional 
pairing (see also the review by Altland et al~), as well 
as by Guruswamy et al 43 for a model with chiral dis- 
order (C z in our notation). Very recently, Aleiner and 
Efetov2& returned to such a RG in the context of disor- 
dered graphene. The renormalization of the conductivity 
gives rise to its dependence on energy (or, equivalently, 
on the electronic concentration, see below). 

Let us briefly analyze the leading logarithmic correc- 
tions and the RG results and compare them to the SCBA 
in the simplest case of diagonal disorder with the only 
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parameter a = 2at r = «o [see Table H] . The diagrams 
of Fig. 01 give logarithmic corrections proportional to 
a log(A/|e|) and missed by SCBA 



5a = 



2e 2 
7r 2 a 



f+2a log(A/| £ |), (a); 
1-2aolog(A/|e|), (b); 



(39) 



(a) 



(b) 



(c) 



(d) 



FIG. 5: One-loop RG diagrams responsible for the renormal- 
ization of (a) the energy and (b, c, d) the disorder couplings. 



A contribution from the diagram (c), which is potentially 
of the same order, vanishes after the angular integration. 
Since the two contributions in H39() cancel each other, the 
SCBA turns out to give the leading logarithmic correc- 
tion even with a correct numerical coefficient, 



a(s) = 



2e 2 
ir 2 ao 



A 

1 - 2a log — 

£ 



(40) 



This coincidence in the numerical coefficient seems to be 
accidental, however. If one takes into account disorder 
amplitudes other than ao, the numerical coefficient in 
front of the leading logarithmic correction in SCBA be- 
comes in general different from the correct one (given by 
RG). 

With lowering energy, consideration of the first-order 
logarithmic correction l|4UI) becomes insufficient. As have 
been already mentioned, all logarithmic corrections to 
the density of states and conductivit y c an be summed 
up with the help of the RG, Refs. [30M4M59ll&cll6lTl6^l6^ . 
Below we briefly present the results in the simplest case of 
long-range Gaussian disorder. For short-range disorder, 
the consideration is similar but five running couplings 
(first 5 in Table [IJ characterizing the disorder should be 
taken into account. As found in Ref. y(J, this does not 
affect qualitatively the behavior of the conductivity. 

After the disorder averaging, the action for electrons 
in graphene with long-range disorder reads 



d 2 r 



(41) 

The vector superfield ip contains 4 x 2 x 2 = 16 compo- 
nents: the four-dimensional structure of the one-particle 
Hamiltonian is complemented by the advanced-retarded 
(AR) and the supersymmetric (boson-fermion) struc- 
tures. The latter serves to perform the disorder averag- 
ing; alternatively, one can use the replica trick. Further, 
A is the third Pauli matrix in the AR space, and the con- 
jugated field is ip = tjj + A. Under renormalization, the 
energy and the disorder strength become running cou- 
plings, e = e(L) and ao = ao(i), where L is the running 
ultraviolet cut-off length (measured in units of vq/A). As 
usual, after elimination of large momenta, the real-space 
coordinates are rescaled to maintain the ultraviolet cut- 
off Vo/A. The coefficient vq of the kinetic term is kept 
fixed by the field renormalization (absent in the one-loop 
order considered below). The relevant one-loop diagrams 
are shown in Fig. [S] (the first two diagrams renormaliz- 
ing ao cancel exactly if the disorder is long-range); the 



resulting RG equations read 
dao 



dlogL 

di 
dlogL 



= 2a 2 , 



= (1 + &o)e. 



(42) 
(43) 



Note that these equations are different from those for the 
random mass problem, Refs. IHlfifiFl only by a sign m 
Eq. f4"2"Jl. The RG equation l|4^|) for the random scalar 
potential problem can be found, e.g., in Ref. l6ll 
Solving these differential equations, we get 



on = 



1 — 2ao log L ' 
sL 

v 7 ! - 2a Q logi 



(44) 
(45) 



The renormalization proceeds until the renormalized en- 
ergy e reaches the value of the cut-off A. Using this 
condition, we find the value of L at which the RG stops, 



A A 
L= Jl-2a log R . 



(46) 



pe. 



The density of states p scales as e L , i.e. peL~ 
Thus, according to Eq. (|45|l . its renormalized value p is 



P 



pL\J\ — 2ao log£. 



(47) 



At the end point of the RG we have p = A/ttvq. Extract- 
ing p from Eq. I|47|l and substituting L from Eq. (|46|l . we 
find the density of states 



Pie) = A 



TTVn 



A 

1 - 2a log — 



(48) 



Further, the conductivity is determined^ by the renor- 
malized dimensionless strength of the disorder ao, 



a(e) 



2e 2 



2e 2 
ir 2 ao 



I - 2a log ■ 



A 



(49) 



in agreement with Ref. l30t 

We see that the result of the SCBA, Sec. lIIIBll agrees 
qualitatively with the fully controllable (RG) solution: 
the conductivity decreases logarithmically up to an expo- 
nentially small scale. The SCBA fails, however, to give a 



8 



correct numerical coefficient in the exponent of Eq. (1331) : 
the correct low-energy scale T is 



T = Ae~ 1/2a °. 



(50) 



Below this new energy scale, the density of states satu- 
rates at a finite value, and the Drude conductivity (with 
localization effects discarded) is of the order of e 2 /h. 
Both these important features are correctly reproduced 
by the SCBA. 

In the experiment, one changes the chemical potential 
/j by varying the gate voltage V g . The electron concen- 
tration n e is proportional to V g , en e — CV g , where C is 
the corresponding capacitance per unit area. Therefore, 
the experimentally measured dependence cr(V g ) is essen- 
tially <j(n e ), up to a simple rescaling. To compare the 
theory with the experiment, we find the density 



n e (ib) = 2 [ de p(e) ~ — : — r. — r. 

VW Jo ™o l-2a logA/| M | 

Combining this with 1)49(1 . we get 



a(n e ) 



2e 2 
ir 2 ao 



1 - a log ■ 



(51) 



(52) 



We see that the dependence of conductivity on electron 
density is only logarithmic, which should be contrasted 
with a much stronger, approximately linear, dependence 
observed in the experiments?^ As we will see in Sec. 
IIII CI such a strong dependence does arise theoretically 
in the limit of strong scatterers. 

One can use a more general RG approach in the case 
of generic Gaussian disorder when all 9 parameters from 
Table |U are non-zero. In doing so, one has to calculate 
the diagrams from Fig. with all possible matrices at 
the vertices of impurity lines. The full set of one- loop 
perturbative RG equations can be found in Appendix IA1 

We note that the conductivity calculated above is not 
the total conductivity far from the degeneracy point. 
There are also weak-localization corrections to the con- 
ductivity, which are small for strong enough dephasing 
or for small enough systems. As shown in Ref. 1261 it is 
convenient to decompose the retarded-advanced Cooper- 
ons in the singlet/triplet representation in both E and A 
channels. Then only singlets with respect to S matter. 
The general expression for the weak-localization correc- 
tion valid for arbitrary disorder then reads 



-^log(^f) [CO 



2c 



(53) 



where / is the electron mean free path determined by the 
renormalized disorder and density of states, and Lir is 
the infrared cutoff set by either the system size or the de- 
phasing length. In Eq. (|53|l . one has to put c; = 1 if dis- 
order preserves the TR-invariance Tj and Cj = (mean- 
ing that the corresponding Cooperon modes are gapful) 
otherwise (see Table QJ. For a combination of several dis- 
order types, only those Cooperon modes remain gapless 



that correspond to the TR-symmetries preserved simul- 
taneously by all disorder matrices involved. In particu- 
lar, for the diagonal disorder ceo ah <k = 1 which yields 
antilocalization, whereas, e.g., for the combination of 7j_ 
and j z disorders we have Co = c 2 = I and c±_ = leading 
to the absence of the one-loop correction. On the other 
hand, for the combination of, e.g., f3 z and "f z disorders, 
Co = 1, c z = c± = 0, and we get localization. Note that 
the weak-localization correction is universal and depends 
only on the symmetry of the HamiltonianS 



C. Unitary limit 



In Sec. IIII Bl we have analyzed the behavior of the 
density of states and the conductivity in the case when 
impurities are weak, so that the disorder can be con- 
sidered as Gaussian. (In terms of the action, Eq. i|41|l . 
this amounted to keeping, after the ensemble averaging, 
only the {^i)) 2 term and neglecting all higher-order cou- 
plings). In this subsection, we will consider the opposite 
case, when the electron is strongly scattered by an im- 
purity and one has to deal with the complete T-matrix 
(|14fl . The analysis of the location of the "phase bound- 
ary" between the domains of weak and strong scatterers 
in the space of microscopic parameters of the problem is 
postponed to Sec. IIII PI 

We proceed by first analyzing the results in the frame- 
work of the SCTMA, Sec. II III and then discuss its ac- 
curacy and limitations. Like in the weak-scatterer limit, 
the SCTMA can be simplified in the limit of strong scat- 
terers. Specifically, at large U we neglect unity in com- 
parison with Ug(s) in the denominator of Eq. (|19|) and 
obtain the following self-consistency equation 



S( £ ) 



77A 2 



[e -£(£)] log_ 



— A 2 



(54) 



The parameter r\ is the dimensionless concentration of 
impurities defined as 



''imp U Q 



1, short-range; 
4, long-range. 



(55) 



The scattering amplitude U does not enter Eq. I|54|) . This 
means the impurities arc effectively considered infinitely 
strong; the limit that we will term the self-consistent uni- 
tary approximation (SCUA). For this type of impurities, 
the weak disorder assumption means that their concen- 
tration is small, rj <C 1. 

The characteristic energy scale in the unitary limit is 
set by the value of £ at zero energy: £(e = 0) = T^ri, 
which we find to be 



A 



V 



l0g(l/7j) ' 



(56) 



In contrast to its Born-limit counterpart T, which is ex- 
ponentially small for a« 1 , the energy scale depends 
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on the disorder strength rj in the power-law fashion. As 
we see below, this is intimately connected with a qualita- 
tively different dependence of conductivity on the Fermi 
energy. 

The further analysis of Eq. (|54|l can be performed in 
the way analogous to our treatment of the Born limit, 
Eq. IMl'l) . We get (see Ref. |58| concerning the crossover 
between high- and low-energy regimes) 



E(e) 



77A 2 



2e 



2(^2 - r») 



e< r 



ii- 



nr sgn e 



log(A/| £ |) ^ 21og 2 (A/| £ |) 



£> r. 



(57) 

Here upper (lower) signs correspond to the retarded (ad- 
vanced) self energy. 

Using Eq. H57|) and the relation between the Green 
function and the self energy in the unitary limit, g = 
— rjA 2 /4nv 2 T,, we get for the density of states, Eq. ipDjl. 



TjA 2 



PSCUA(e) = — 5— 2I ImS 1 (s)\ - ; 



e > r J( . 



The density of states is constant at small energy and 
shows the linear dependence characteristic for the clean 
graphene at high energies. To find the disorder correc- 
tion to this result, one has to use a more precise value of 
the self energy than that given by Eq. i|57|) . After two 
iterations of the equation (|54|l . the linear- in- 77 contribu- 
tion to the density of states is obtained (see details in 
Appendix lB|l . 



(1) / \ 

PSCVA\ £ ) 



[1-aiKe)] 



(59) 



Here the parameter ajj has the meaning of the inverse 
dimensionless conductance [see Eq. (|63ll below], 



au{e) 



77A 2 



''imp 



A 2 



2 £ 2 log 2 (A/| £ |) log 2 (A/| £ |) 



(60) 



It is of the order of the squared ratio of the electron wave- 
length A e at energy s to the distance between impurities, 
up to a logarithmic factor. The condition e 3> T n ensures 
that the relative correction is small. 

To find the correction to the density of states of the 
second order in 77, one has to go beyond the self-consistent 
approximation. The diagrams with intersecting impurity 
lines become important in this case as we have already 
seen it in the Born limit (Sec. IIII B 2() . A rigorous cal- 
culation taking into account all second-order diagrams is 
given in Appendix [5] The result is 



A 



A 



1 - au(e) - 2a(j(e) log — log log — 



(61) 



In order to calculate the conductivity, we find the po- 
larization operator, Eq. (|23[) . 



U RA {e) = 



1 r)A 2 £ -2RcS( £ ) 
4™o 2|E( £ )| 2 e-Re£( £ ) ' 



(62) 



Substituting it into Eq. I|25|l . we obtain for the conduc- 
tivity at not too low energy, e T v , 



o'scuA(e) — 



2_2 



4e 2 £ 
TT 2 r)A' 



, 2 A 



(63) 



Equation l|63|) is written for the case of long-range disor- 
der; if the disorder is short-range, the vertex correction 
is absent and the resulting conductivity is twice smaller. 

What about the multiple scattering of electrons on 
complexes of two or more impurities (described by mutu- 
ally "entangled" T-matrices, Appendix [0) that are not 
included in the SCUA? We remind the reader that in 
the Born limit of weak impurities fSec. IIII B 2() . similar 
multiple scattering processes contribute to the dominant 
(logarithmic) energy dependence of the conductivity, Eq. 
(|49|l . In the unitary limit, however, the dominant energy 
dependence of the conductivity a(e) cx l/au(e) comes al- 
ready from the £ -dependence of a T-matrix describing the 
scattering off a single impurity. Therefore, the logarith- 
mic corrections to the conductivity, analogous to those 
in the Born limit [see Eq. (|49l) ]. are of minor importance 
in the unitary limit. 

As in the case of Born-type disorder, Sec. 



IHTBI we 

now convert the energy dependence of conductivity into 
its dependence on the electron concentration n e . We have 
according to Eq. (for /j » r ?) ) 



n e (/i) 

so that Eq. yields 
a(n e ) 



At: 2 n- u 



irv 2 



log 



A 2 



v 2 \n e \' 



(64) 



(65) 



For the short-range disorder, the result is twice larger. 
In contrast to the limit of weak (Born) scatterers, the 
conductivity shows a strong concentration dependence: 
it varies linearly with n e , with a logarithmic correction. 
This result compares nicely with the experimentally ob- 
tained linear behavior of cr(n e ) (or, equivalently, constant 
mobility) m& This indicates that the dominant scatterers 
are strong. Equation Ij65(l predicts a logarithmic correc- 
tion to the linear behavior, which should become more 
pronounced if the measurement is extended to larger gate 
voltages. 

One can also calculate the SCUA conductivity at 
e <C T,,. In this limit, the contributions (jG R jG R ) and 
(jG A jG A ) should also be taken into account. The Drudc 
conductivity then appears to have exactly the same value 
ctscua = 2e 2 /n 2 h as in the SCBA, Sec. IIII B II Analo- 
gously to the SCBA case, this result is questionable in 



10 



view of the localization effects neglected in the Drude 
formalism. It is important to recall in this context that 
infinitely strong impurities are chiral (Cr), yielding a di- 
vergent density of states^ (cf. Refs. l50ll5lT) in this situ- 
ation. We will return to the conductivity at half filling 
for a chiral disorder in Sec. IIVI 



D. Phase diagram 



In the preceding subsections, Sec. IIII Bl and IHlUl we 

have studied the limits of weak (Born) and strong (uni- 
tary) scatterers, respectively. We have found that the 
behavior of the conductivity is essentially different in the 
both limits: it depends only logarithmically on energy 
in the Born limit, and shows a linear behavior (with a 
logarithmic correction) in the unitary limit. The aim of 
the present subsection is to construct a "phase diagram" 
that would predict which of these types of behavior is ex- 
pected for given characteristics of disorder. (Of course, 
we do not mean any phases in the strict sense; there is 
a smooth crossover between the Born and the unitary 
regime) . 

The unitary limit corresponds to the neglect of unity 
in comparison with Ug(e) in the denominator of Eq. (|19J) . 
In order to see when this is justified, we use the large- 
energy expression i|57|) for T,(e) and compare Ug{e) with 
1. This yields the following energy-dependent value of 
the parameter r\ at the "phase boundary" between the 
weak-scatterer and strong-scatterer regimes, 



i \ a£2 i 2 A 



(66) 



or, equivalently, olu(e) ~ a. For r\ <C »7c( £ )i that is 
au(e) <C a, the system is in the unitary limit. 

In Fig. |21 we plot the phase diagram in both e-rj and 
e— a coordinates. Remarkably, the system may pass from 
the Born into the unitary limit when the energy in- 
creases while the disorder remains fixed. In Fig. this 
crossover occurs for a broad range of impurity concentra- 
tions, a~ 1 exp(— 1/a) <C i] « a. At still smaller values 
of r\ the system is in the unitary phase at all energies. It 
is worth stressing that the unitary phase in Fig. HJ is es- 
tablished even for a -C 1, when disorder could be naively 
considered as weak. The reason for this effect is as fol- 
lows. The growth of the density of states with increasing 
energy results in a more efficient scattering of higher en- 
ergy electrons by an impurity, thus making the scatterer 
effectively stronger at higher energies. With increasing 
a the phase boundary l|66|) in Fig. Amoves upwards, and 
for a > 1 the Born phase disappears altogether. 

A unitary-to-Born crossover discussed above would 
manifest itself in a change of the behavior of the con- 
ductivity, from the linear energy dependence at high e 
(Sec. IIII Bl) to the logarithmic dependence at lower en- 
ergies (Sec. nuTty . 66 Experimentally, the measured con- 
ductivity of graphene shows a linear dependence down 
to the lowest energy scale (where a saturates at a value 



a 



Born 



unitary 



o r 





unitary 






Born 



FIG. 6: "Phase diagram" in the £-77 plane for a fixed a (upper 
panel) and in the e— a plane for a fixed 77 (lower panel). The 
solid line is the "phase boundary", Eq. 16611 or Eq. I6UII . where 
the crossover between the Born and the unitary regimes takes 
place. At low energies (dashed part) the density of states 
saturates, while the Drude conductivity reaches a value ~ 
e 2 /h, implying that generically the localization effects should 
become strong. 



~ 4e 2 //i). This indicates that the scattering is domi- 
nated by strong impurities, which remain in the unitary 
part of the phase diagram down to the lowest energies. 



E. Charged impurities 

The case of charged impurities deserves a special con- 
sideration. If such impurities are located far from the 
graphene layer, they are expected to be screened by the 
gate and will not be different from finite-range scatter- 
ers considered above. Let us consider, however, charged 
impurities located near the graphene layer. The scatter- 
ing potential of the Coulomb center in 2D is Vb(q) = 
27re 2 /x<Z- Taking into account the static screening by the 
graphene electron gas in the random phase approxima- 
tion (RPA), we obtain^ 



X q + 2ire 2 p(e) ' 



(67) 
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where \ is the dielectric constant. Strictly speaking, the 
RPA is not justified in graphene since the parameter r s — 
e 2 /Hvox is of order unity. It will be sufficient, however, to 
find a parametric behavior of quantities under interest, 
up to numerical coefficients of order unity. 

As follows from Eq. (|67|l . the intervalley-scattering 
component of the Coulomb potential, V(ko) ~ 2ne 2 /xko 
is very small and can be neglected, so that the Coulomb 
impurities are of long-range type. As to the scatter- 
ing within one valley, it is only slightly anisotropic. In- 
deed, the inverse screening length k = lite 2 p{e) is of the 
same order that the characteristic momentum transfer, 
k ~ q ~ e/vo for r s ~ 1. Therefore, up to a numerical 
factor of order unity, we can neglect q in the denomina- 
tor of Eq. i|67|) (which means a neglect of the anisotropy 
of the intra- valley scattering) . This brings the screened 
charged impurities into the class of long-range scatterers 
considered above but with an energy-dependent ampli- 
tude, 

9 

U{e) = p-\e)^^-. (68) 

There is, however, an important difference between a 
charged impurity and a long-range potential impurity. 
The scattering amplitude for slow electrons, with mo- 
menta q < |e|/vo is given by Eq. (|rT8l . while the electrons 
with larger momenta arc scattered much less efficiently 
due to the lack of screening at small distances. This can 
be taken into account by setting an effective high-energy 
cut-off A ~ |e|. 

Finally, using Eqs. (|3Uj> . (|55j> . and l|6l)jl with U from Eq. 
(I68JI and A ~ e, we obtain -q c ~ r\. Thus we come to the 
conclusion that, with charged impurities, the system is 
just at the crossover between Born and unitary regimes. 
This also justifies the use of the clean density of states in 
Eq. I|68|) . Indeed, approaching the crossover from unitary 
side, the disorder-induced corrections to p(e) are negligi- 
ble, see Eq. (|59|l . On the other hand, if one uses the Born 
expression Eq. (|48|l . logarithmic corrections are absent, 
as long as A ~ e, and the clean value of the density of 
states in Eq. (|B^|l is again justified. 

The energy and density dependencies of the conductiv- 
ity of graphene with Coulomb impurities are thus equiv- 
alently given by both Born [Eqs. (|49|l . (|52|) with energy- 
dependent coupling ao ~ ni mp v 2 /e 2 ] and unitary [Eqs. 
(JHSl: (E3] expressions with logarithms omitted, 

(69) 

^imp^o ^imp 

The Born approximation was use d for c alculating the 
conductivity in recent works Refs. Il6l28l (see also Ref. 
13^) . The result is consistent with Eq. I|69|l . A different 
result (containing an additional logarithmic factor) was 
obtained in Ref. [7CJ. We believe that the derivation in 
Ref. is incorrect since it employs the quasiclassical 
Thomas-Fermi approximation beyond its range of valid- 
ity (at energies much larger than ep). 



IV. CONDUCTIVITY AT THE DEGENERACY 
POINT: CHIRAL DISORDER 

A. Universal conductivity 

In this section we consider the conductivity of graphene 
at half filling, e = 0. The Drude conductivity obtained 
self-consistently in Sec. lIIIl in both Born and unitary limit 
has the value a = 2e 2 /ir 2 h at this point. Since this value 
is of the order of conductance quantum, this is by no 
means the end of the story: the localization effects be- 
come strong at half filling. If the intervalley scattering 
is weak (long-range disorder potential), an intermediate 
temperature range exists where the conductivity correc- 
tion is positiva2ii2§i22iZ£ due to the additional Berry phase 
7r associated with the electron pseudospin in the sublat- 
tice space. This situation belongs to the symplectic sym- 
metry class. With lowering temperature T, the interval- 
ley scattering comes into play and a crossover to the or- 
thogonal symmetry class occursi£§i22*2I The localization 
correction becomes negative and drives the system into 
the strong localization regime^ Thus, for a generic dis- 
order, the conductivity at half filling should have a pro- 
nounced temperature dependence and get strongly sup- 
pressed with lowering T. Surprisingly, this is not what is 
observed in the experiment. The conductivity has been 
founds to be close to the value 4e 2 /h, remaining T- 
independent in a broad range of temperatures. The aim 
of this section is to analyze whether and in what situation 
this behavior may be expected theoretically. According 
to what was said above, this might only happen, if at all, 
for a particular type of disorder. 

The special class of disorder that we will consider in 
this section is the randomness that preserves one of the 
chiral symmetries <|13f) of the clean graphene Hamilto- 
nian. Some possible realizations of such type of disorder 
were listed in Sec. Ill CI Whether the dominant disorder 
in graphene may be of this kind is an open issue, which 
may be related to technological aspects of the sample 
preparation. Our aim here will be to analyze what are 
consequences of the assumption of chiral character of dis- 
order. 

A peculiar behavior of 2D systems with chiral disor- 
der with respect to localization effects has been demon- 
strated by Cade and Wegneri 67 ! 73 They considered a ran- 
dom hopping problem on a square lattice and showed 
that at zero energy, where the system possesses the chi- 
ral symmetry, the RG /3-function of the corresponding 
cr-model vanishes to all orders in the inverse conductiv- 
ity, implying that the conductivity is not renormalized. 
This absence of usual infrared-singular corrections to the 
conductivity due to cooperon- and diffuson- loops can be 
attributed to the fact that the "antilocalizing" interfer- 
ence corrections to the density of states cancel the lo- 
calization corrections to the diffusion coefficient. The 
density of states has been found^SLZ^i 7 ^ to diverge as 
p(e) ~ e~ 1 f(e) for e — > 0, where /(e) gives the sublead- 
ing e-dependence and provides the convergence of the to- 
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tal number of electronic states^ At any finite e the chi- 
ral symmetry is broken and the localization on the scale 
£(e) oc |/(e)| _1 ^ 2 occurs^ The states at the band cen- 
ter e = are delocalized and the conductivity a[e = 0) 
takes a finite value depending on the disorder strength. 
According to the cl assifica tion of Refs. l63l76l77t the sys- 
tem studied in Refs. I67l73l belongs to the chiral symmetry 
class AIII. 



While the results of Refs. 67 73 suggest that one may 
expect a finite zero-energy conductivity in our problem, 
they cannot be directly applied. Indeed, the dimension- 
less Drude conductivity at e = is of order unity in our 
case, whereas it should be large to justify the derivation 
of the a model and of the perturbative RG. Another re- 
lated peculiarity of the problem we are considering is the 
Dirac dispersion of carriers. This will allow us to prove 
below a statement that is still stronger than that of Gade 
and Wegner: we will show that for certain types of chiral 
disorder all disorder-induced contributions to conductiv- 
ity cancel. 



1. Co-chirality: symmetry consideration 

Let us consider the disorder which preserves the Co- 
chirality, H = — (i^Has. The random part of the Hamil- 
tonian contains matrices 0OT3, Ci^ti^, and ui^tq- Ac- 
cording to the Table I, in the case of weak disorder, the 
corresponding coupling constants are a±, j3±, and jx- 
While the disorder characterized by (3±_ and 7j_ preserves 
the time-reversal invariance To, the a^-disorder, being 
physically a random vector potential, violates the To _ 
symmetry. 

According to Ref. |7^, the random Dirac Hamiltoni- 
ans preserving the Co-chirality and violating the TR- 
symmetry (case 1 of Ref. l78f) belong to the chiral symme- 
try class AIII, while the combination of Co-chirality and 
To-symmetry (case 6 of Ref. drives the system into 
the Bogolyubov - de Gennes symmetry class CI. In both 
cases, the low-energy theory (a model) is affected by the 
presence of the Wess-Zumino-Witten term in the action. 

The one-loop RG equations for Co-disorder read (see 
Appendix E}: 



da± 
dlogL 

d(3± 
dlogL 

<91ogL 



0. 

4/3_L7±> 
Pi- 



(70) 
(71) 
(72) 



Note that the equation (|70|l for a± is split from Eqs. i|7ip. 
(I72|l . This set of equations is identical to Eq. (19) of Ref. 
|63j with g' = aj_, g = 2j±, g v0 = 0, and g^ = 2(3±. In 
Ref. 163, these couplings described the scattering between 
the four nodal points of the spectrum of disordered d- 
wave superconductor. Our problem with only two nodes 



corresponds to setting the coupling between the neigh- 
boring nodes in d-wave superconductor to zero, g^o = 0, 
while retaining the forward-scattering (intranode, g) and 
backscattering (scattering between the opposite nodes, 
g-K-x) amplitudes. This situation (non-Abelian vector po- 
tential problem) was considered in Ref.0 (see also Refs. 
iTlillSGT) . where the density of states was shown to vanish 
in the limit e — > as 



p(e) ex lei 1 / 7 . 



(73) 



Here 1/7 = 1/(2N 2 - 1), where = 2 is the number of 
flavors (nodes). 

In the presence of the random vector potential only 
(aj_ coupling, preserving all the four chiralities simulta- 
neously, class AIII), the density of states also goes to 
zero with decreasing energy, but with a non-universal 
exponent 60 61 which depends on a±: 



p(e) cx | £ |(1-«l)/(i+«l). 



(74) 



Note that in this random vector potential problem, the 
disorder strength remains non-renormalized, see Eq. I|7U|) 
[in fact, the one-loop equations for a± and e are exact, 
see Refs. 43.60]. Therefore, a± does not generate the 
scale r and the one-loop result for the density of states, 
Eq. 174|) . holds in the whole range of energies below A. 



2. Co-chirality: conductivity at Dirac point 

We are now going to study the conductivity in the 
situation when disorder preserves Co-chirality. The chiral 
symmetry Co allows one to relate retarded and advanced 
Green functions: 

<7 3 G*< A >( e ;r,r>3 = -G A ^\-e;v,v'). (75) 

The conductivity is given by the Kubo formula l|21|l . 
which we rewrite here in the full form, 



1 



d 2 (r - r') Tr j x G R (0, r, r')j x G A (0, r', r) 



--fG fi (fl,r,r')fG fi (0,r',r) 

-~fG A (0,r,r')fG A (0,r',r) 



(76) 



Now we use the identity l|75(l to trade all advanced Green 
functions in Eq. i|76|) for retarded ones and thus to present 
the conductivity in terms of retarded Green functions 
only. Further, we exploit the following important relation 
between the components of the current operator 



&3j 



-foz = if, 



(77) 



which is the consequence of the Dirac spectrum. At this 
point, our problem differs from that considered by Gade 
and Weg 

ner 67,73 

who dealt with a bipartite square lattice 
with a non-linear electronic spectrum. 
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The transformations Eqs. (|75(l and l|77(l allow us to cast 
the Kubo formula in the following form: 

°™ = -\ E Jd\r-r>) 

■[j a G R (0;r,r')j a G R (0;r',r)\. (78) 



a—x,y 

x Trl 



At first glance, this expression is zero due to the gauge 
invariance. Indeed, the right-hand side of Eq. I|78|) is pro- 
portional to the second derivative of the partition func- 
tion Z[A] = TrlogG^A] (or, equivalently, first deriva- 
tive of the current Tr j a G R [A]) with respect to the con- 
stant vector potential A. The gauge invariance implies 
that a constant vector potential does not affect gauge- 
invariant quantities like the partition function or the 
current, so that the derivative is zero. This argument 
is, however, not fully correct, in view of a quantum 
anomaly present in this problem. The elimination of A 
amounts technically to a shift in the momentum space 
k — * k — eA, which naively does not change the momen- 
tum integral. If we consider a formal expansion in the 
disorder strength, this argument will indeed hold for all 
terms involving disorder but not for the zero-order con- 
tribution. The momentum integral J cPk Tr j a G R (k) is 
ultraviolet-divergent and the shift of variable is illegiti- 
mate. This anomaly was first identified by Schwinger 81 
for 1 + 1-dimensional massless Dirac fermions. In the 
Schwinger model, the polarization operator is not af- 
fected by an arbitrary external vector potential A(ir, t) 
and is given by the anomalous contribution, yielding a 
photon mass in the 1 + 1 electrodynamics^^ In our 
analysis, the role of A(x, i) is played by the chiral disor- 
der. The explicit calculation of the zero-order diagram 
(the one with no disorder included) yields 



ie 2 v 2 



d 2 k 



ir J (2tt) 2 {v 2 k 2 + S 2 ) 2 



2e^ 

IT 2 ' 



(79) 



Here 5 is an infinitesimal imaginary part in the denomi- 
nator of the Green function; we will return to its role and 
physical meaning below. We note that the same universal 
value of the conductivity in the situation when the only 
type of disorder is the abelian random vector potential 
(a±_) was previously obtained in Ref. l60l 

An alternative derivation of the same result is based 
on the Ward identity 



- le (r-r')G fl (0;r,r') = [G R iG R ](0;r,r'). 



(80) 



Averaging it over disorder, plugging it into Eq. (|78|) . 
transforming to the momentum space, and performing 
the integration by parts, we are left with the surface con- 
tribution only, 



* = -^fdk n Tr[jG R (k)]. 



(81) 



where the integral is taken over a large circle |k| = 
const — > oo. For large momenta the Green function can 



be replaced by its bare value, which yields again the uni- 
versal conductivity 



£_ I dknk - 

~3 f £-2 ~ 



7T 



(82) 



This universal value of the conductivity is independent 
of the ultraviolet cut-off in the momentum space. This 
signifies that the integral in Eq. I|78|l is accumulated in 
the vicinity of the degeneracy point, as seen explicitly 
in Eq. (J73J). The fact that, in a realistic system, the 
linearization of the spectrum ceases to be valid at high 
momenta does not spoil the derivation: the functions G R 
and G A are essentially equal to each other there, so that 
the integrand of Eq. (|76|1 is cancelled. 

It is worth emphasizing that, as is clear from the 
derivation of Eq. H82|) . it assumes that the ultraviolet 
cut-off A is much larger than the disorder-induced en- 
ergy scale r. (More accurately, here T is the low-energy 
electron relaxation rate determined as a scale where the 
dimensionless Drude conductivity is of order of unity, or, 
equivalently, where the RG flow enters the strong cou- 
pling regime.) In other words, the disorder is weak, i.e. 
a C 1 for Gaussian disorder. One more formulation of 
this condition is that for energies comparable to the cut- 
off, e ~ A, the Drude conductivity considered in Sec. 
IIIII is large (compared to e 2 /h). This condition, that we 
assume throughout the paper, is very well fulfilled in the 
experiments^*^ Violation of this condition would imply 
that the disorder is so strong that it completely destroys 
the Dirac character of the spectrum. In this situation the 
universal value of the conductivity (|79H , H82f) of the chiral- 
symmetric system would not survive. The corrections to 
the universal value of the conductivity are exponentially 
small: 8a < e 2 T/A, which implies that there are no cor- 
rections to any order in the perturbative expansion of 
a(e = Q) in a < 1. 

The above derivation of the universal conductivity re- 
mains valid for the case when a magnetic field of an arbi- 
trary strength is applied: the vector potential A a couples 
to the current, i.e. to the matrices a a , a = x,y, thus pre- 
serving the chiral symmetry. In this context, it is worth 
mentioning the result of Hikami, Shirai, and Wegner— 
who found that the longitudinal conductance in the cen- 
ter of the lowest Landau level of the chiral-disordered 2D 
electron gas is equal exactly by a — 2e 2 /ir 2 H in the limit 
of very strong magnetic field, when the Landau level mix- 
ing can be neglected. Their finding can be considered as 
a B — > oo limit of our general result. Indeed, in this limit 
the kinetic energy is frozen, so that the difference between 
the electron dispersion on the square lattice (considered 
in Ref. l83|) and the graphene lattice becomes immaterial. 

We turn now to an important and delicate point related 
to the above derivation of the universal conductivity 17911 . 
Specifically, we have introduced an infinitesimally small 
imaginary part of energy, 5. Physically, it has a mean- 
ing of the electron lifetime or, alternatively, a dephasing 
rate, and can be thought as modelling processes of es- 
cape of electrons in some reservoir or some dephasing 
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mechanism. Models with such a uniform constant value 
of 5 were used in the literature to imitate dephasing in 
quantum dots, see e.g. Ref. 0. 

In our calculation, S has served as an infrared regulator 
for the theory. Although it has dropped from the final 
result, its role is not completely innocent. Depending on 
the physical situation, the infrared regularization may be 
provided by different quantities, which, as we are going 
to discuss, will influence the value of the conductivity. 
Specifically, in addition to 5, we can imagine the follow- 
ing sources of the infrared cut-off: (i) finite frequency, (ii) 
finite system size, and (iii) interaction-induced dephasing 
at finite temperature. In Sec. IIV Bl we will analyze the 
frequency dependence of the conductivity. As to the sit- 
uations when the temperature or the system size govern 
the infrared behavior, we restrict ourselves to brief com- 
ments only, relegating a detailed analysis to future work. 



3. C z -chirality 

Let us now turn to the disorder which preserves the 
C z -chirality, H = —(73X3^(73X3. The random part of the 
Hamiltonian may then contain matrices 173x1,2, 01,2X3, 
0i, 2X0, and <jqT\^,- The first two (the corresponding cou- 
pling constants are 0q and a±_) violate the time-reversal 
symmetry To, the last two (7j_ and Z ) preserve it, see 
Table I. Note that the disorder characterized by a± and 
7_i_ (real and imaginary vector potentials, respectively) 
also preserves the chiral symmetry Co considered above. 

According to Ref. l78l random Dirac Hamiltonians pre- 
serving the C^- chirality and violating the TR-symmetry 
(case 2 of Ref. l78j) belong to the chiral unitary symmetry 
class AIII. The combination of C z - chirality and the time 
reversal invariance To (case 9+ of Ref. l7ql corresponds 
the chiral orthogonal symmetry class BDI. Finally, the 
combination of C z -chirality and T z -symmetry (case 9_ of 
Ref. ITsI) falls into the chiral symplectic symmetry class 
CII. 

The one-loop RG equations for C z -disorder read (see 
Appendix [XJ : 



da± 
dlogL 

dlogL 
O0 Z 

dlogL 
dj± 

dlogL 



This model was considered in Ref. |4J; the RG equations 
iJBBJ-JHU agree with the set of equations (4.84) in Ref.li^ 
with g^ = O /V2, g A > = a±, 9A = 7-U and g m = (3 x /y/2. 
If the system is time reversal (To) invariant, only the 
coupli ngs Z and 7 Z survive; this case was considered in 
Refs. 43 85. The density of states in the generic C z -case 



20oPz, 




(83) 


2a ± {0 o - 


~0z), 


(84) 


2a ± {0 o ~ 


~0z), 


(85) 


0o+0l 




(86) 



divergea 43 i 67 i 68 i 73 i 74 i 75 in the limit e -»■ 0, see Sec. ITVAl 
[for the case of the random vector potential, see Eq. (|71|) ]. 

Let us turn to the conductivity at half filling for a 
generic disorder preserving the C z chirality. The proof 
of the universality of the conductivity based on gauge- 
invariance argument does not work now. Indeed, the 
C z -chirality transformation of the Green's function 



03T3 



e;r,r Jct 3 x 3 



-G A{R \-e;r,r') (87) 



generates the new vector vertices j x ' v Ts instead of cur- 
rents (these new vertices can be considered as J5 currents 
of Dirac fermions). Then we are left with the G R G R - 
type correlators of both j x ' v and j x ' v ■ The latter can 
not be obtained as derivatives of the partition function 
with respect to the constant vector potential A. Never- 
theless, for weak disorder we find that the conductivity 
at half filling is still universal, a(e = 0) = 2e 2 /7r 2 , up to 
corrections in powers of disorder strength. 

To show this, we first calculate the perturbative cor- 
rection 5a^ to the conductivity of a pure system at the 
first order in disorder strength and find that it vanishes, 
Sa^- 1 ' = 0. This implies that the conductivity at e = 
does not depend on the ultraviolet cutoff A. Indeed, all 
the contributions generated by the RG (and thus depend- 
ing on the ratio A/S) sum up to zero, because we can use 
the fully renormalized disorder as an effective single im- 
purity line in 8<j^ — 0. The second-order perturbative 
calculation yields 



5a 



(2) 



2tt 2 



(A> - Pzf 



We note that the combination 0o — Z is not renormalized 
during the RG procedure, as follows from Eqs. (|84|l and 
(|85|l . This is in agreement with the above RG argument 
for the first-order correction. Thus the conductivity at 
the Dirac point can be presented as a series in the pa- 
rameter 0o — Z . Next, we recall that for C z -chirality, 
the RG /3-function of the Gade-Wegner a model 67 * 7 ^ van- 
ishes to all orders, so that there are no singular quantum- 
interference corrections to a{e — 0) due to the soft modes 
(impurity ladders). This proves that the expansion of 
<j(e = 0) in powers of O — Z converges. Thus for the 
case of weak disorder the conductivity is universal with 
small corrections in powers of the disorder strength (un- 
like in the case of the Co-chirality, where the corrections 
are nonperturbative in the disorder strength.) 



4- C±-chirality 

Finally, let us discuss the case of C^y-chirality (cou- 
plings ax, 7oj and 7 Z ). Each of these chiralities taken 
separately is similar to the C z -chirality. However, in an 
isotropic system considered here, both C x and C y chiral- 
ities are expected to be present simultaneously. This im- 
plies that the disordered Hamiltonian anticommutes with 
both xi and X2 and hence is proportional to X3 . Thus it is 
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split into two equivalent copies. Therefore, the symme- 
try of the problem is governed by the properties of the 
"sub-Hamiltonians" and its chirality is in fact fictitious. 
In particular, the generic case with all a±, 70, and 7z 
present jSi corresponds to the conventional Gaussian uni- 
tary class A (Quantum Hall effect). A single coupling 
"f z corresponds to the symmetry class D (random mass 
problem). In all these cases the system is in a critical 
phase so one can expect a finite conductivity at e = 0. 
For the sake of completeness we present the RG equations 
for the Cj_-chirality: 



da± 
dlogL 

d-fo 
dlogL 

d 7z 
dlogL 



4 7o7z , (89) 
2(«± +7o)(7o +7*), (90) 
2(a ± -7 2 )( 7o + 7z ). (91) 



We are not aware of realistic examples of the disorder pre- 
serving the C^-chirality in the context of the transport 
in disordered graphene. Therefore, we will not consider 



this case in the rest of the paper. 



B. Conductivity at finite frequency 

In this subsection, we analyze the frequency depen- 
dence of the conductivity. For completeness, we also keep 
a small level width 6 introduced above. It was crucial 
for the argument leading to Eq. I|79|l that the system is 
exactly at half filling, e = 0. A non-zero frequency im- 
plies an integration over the energy range of the width lu, 
which breaks the chiral symmetry. When the frequency 
lu is much smaller than S, this effect is however negligi- 
ble, the infrared regularization is provided by <5, and the 
universal result (|79|l survives. In its turn, 6 plays no role 
when it is the frequency that serves as a dom- 

inant infrared cut-off now. In the high-frequency limit, 
lu r, the situation simplifies again: one can neglect the 
effect of disorder altogether and calculate the conductiv- 
ity by a simple Kubo formula with bare Green functions. 
The result for the real part of the conductivity is again 
universal but with a slightly larger valuei22i2£ 



2 r u de 
Re a (lu > T = - / — 

7T Jo OJ 



d 2 k 



■Tr 



3 x ImG«(e-Lu,k)flmGK(e,k) 
de f d 2 k 



= 87re 2 ^ 



e 

T 



(92) 



A very interesting new situation arises in the interme- 
diate regime, { « w « T. Here lu plays a twofold role, 
leading to two competing effects. On one hand, as dis- 
cussed above, the frequency drives the system away from 
the chiral-symmctric point and thus restores localization. 
On the other hand, the frequency cuts off the singular 
localization correction. Which of these effects wins? To 
answer this question, one should compare lu with the level 
spacing in the localization area, A^(e), where e ~ lu. In 
order to find the scaling of with energy, we consider 
a RG transformation that drives the system away from 
the chiral fixed point. The RG stops when the renor- 
malized energy e reaches the macroscopic scale A; on 
such scales the disorder becomes already strong since the 
initial value of e was below V. In this strongly disor- 
dered case, the value of the running ultraviolet cut-off 
length Lvq/A (corresponding to the renormalized elec- 
tron wavelength) determines then the localization length 
£. As discussed in Sec. IIII B 21 the density of states p 
scales as e _1 L 2 . Therefore, 



implying for the level spacing at the length £, 



(«b/A) s 



pe 



(93) 



A £ (e) 



1 



P (e)e(e) 



(94) 



This result is rather general and is only based on the fact 
that the operator governing the flow of the system away 
from criticality couples to the energy in the action. One 
can of course explicitly verify that the results of Ref. 
[67] for the density of states and the localization length 
quoted in Sec. IIV AI satisfy Eq. 

We conclude that the two competing effects of the 
frequency (the localization and the infrared regulariza- 
tion) "make a draw" - both of them are equally impor- 
tant. Therefore, the system turns out to be, roughly 
speaking, half way between the chiral fixed point and 
the conventional symmetry. This results in a new uni- 
versal (frequency-independent) value of the conductiv- 
ity <Ju ~ e 2 /h in the considered regime S <C lu <C T. 
More precisely, this value depends on the type of chiral- 
ity and the symmetry class of the system away from the 
degeneracy point. In particular, the system with generic 
Cq and C z chiral disorder with (without) TR symmetry 
To is driven into the Wigner-Dyson orthogonal (respec- 
tively, unitary) symmetry class by finite energy. On the 
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2e 2 



FIG. 7: Frequency dependence of conductivity in a system 
with chiral disorder. At intermediate frequency, 8 <C OJ <SC T, 
the conductivity acquires some universal value er^ of the order 
of e 2 /h which is not known analytically. This value depends 
on the type of chirality and the symmetry class of the system 
away from the degeneracy point. 



other hand, the system with C z - and T z -invariant disor- 
der (/3o and 7j_) falls into the Gaussian symplectic sym- 
metry class away from e = 0. 

The frequency dependence of the conductivity is 
sketched in Fig. Despite its universality (for a given 
symmetry) , the value cr w most likely cannot be calculated 
analytically, since this would require an exact knowledge 
of the full crossover between the chiral and the normal 
classes. 



Aj(T) is ~ T, so that the result of the competition of 
these two effects depend on the value of Tr^(T). The 
theory of dephasing in the present situation remains to 
be developed. If the dominant mechanism of dephasing 
is the electron-electron interaction, one can expect that 
(like in conventional 2D systems with dimensionless con- 
ductivity replaced by unity) t^ 1 (T) ~ T. If this is indeed 
true, the behavior of the conductivity at T < T will be 
qualitatively analogous to that for the case of finite fre- 
quency, Sec. IIVBI At high temperatures, T > T, the 
T-dependence of the conductivity will essentially repro- 
duce its e dependence for given type of disorder (Born or 
unitary). 

More realistically, one can think about a situation 
when the disorder is predominantly chiral, but the chi- 
ral symmetry is slightly broken^2<£i e.g., by weak po- 
tential disorder on the energy scale T x <C T. Then the 
above consideration allowing one to expect the conduc- 
tivity ~ e 2 /h will be applicable in the intermediate range, 
r x < T < T; at still lower temperatures, T -C T x , the 
chirality-breaking effects will drive the system into the 
strong localization regime. 

It is worth noting that the interaction may lead to 
other effects (in particular, to open the gap in the spec- 
trum and/or to break the chiral symmetry, cf. Refs. 
I25H80T) not included in our consideration. These ques- 
tions also require further study. 



2. Mesoscopic sample 



C. Additional comments 

In Sees. IIV Al and ITY Bl we have analyzed the conduc- 
tivity in the case when the dominant infrared regulariza- 
tion is provided either by the inverse life time 6 or by 
the frequency u. As has been mentioned above, this role 
may be alternatively played by the interaction-induced 
dephasing at finite temperature or by the finite size of 
the system. Leaving a detailed analysis of these prob- 
lems for the future, we only make some comments on 
them i n Sees. IIV C II and IIV C 21 below. Finally, in Sec. 
IIVC 31 we briefly discuss what happens with the problem 
considered when we pass from the 2D geometry to the 
quasi-lD one by rolling the plane into a cylinder. 

1. Temperature dependence 

In the presence of interactions, the temperature T 
plays a twofold role, similarly to the frequency. On one 
hand, it induces an averaging over the energy window 
of the width ~ T, thus breaking the chiral symmetry 
and "switching on" the localization effects. On the other 
hand, the interaction at finite T generates a non-zero 
dephasing rate t ( ^" 1 (T) cutting off the localization cor- 
rections. As we showed in Sec. IIVBI the level spacing 



Let us now consider the situation when all the potential 
infrared regulators 8, T, to are much smaller than the level 
spacing in the sample. In this case, the sample will be 
fully phase-coherent (mesoscopic) and its size will serve 
as an infrared cut-off. In such a mesoscopic situation one 
should in general speak about a conductance, not conduc- 
tivity. Furthermore, the properties of the conductance 
will essentially depend on the sample geometry. Con- 
sider a rectangular sample L x X L y , with current flowing 
along the x axis. For an approximately square sample, 
L x ~ L y , we expect, based on the above results, an aver- 
age conductance of order e 2 jh. Indeed, one can imagine 
taking first 6 larger than the level spacing (so that the 
result of Sec . IIV Al applies) and then decreasing it until it 
reaches the level spacing.) The above statement follows 
from the continuity. In view of the mesoscopic character 
of the sample, we also expect in this case a broad distri- 
bution of the conductance, with a variance of the order 
of (e 2 /h) 2 . Both the average value and the conductance 
distribution will depend on the exact value of the aspect 
ratio L x /L y in a non-trivial way and can hardly be cal- 
culated. 

For a long sample, L x 3> L y , the geometry becomes 
quasi-one-dimensional, and our results for the universal 
conductivity ~ e 2 jh cease to be relevant (see also Sec. 

Eg. 

Finally, let us consider a case of a very broad and short 
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sample, L y ^>L x . In this situation, the conductance will 
be self-averaging, so that one can again speak about con- 
ductivity. Using again the continuity, we conclude that 
the conductivity in this situation will have some univer- 
sal value ox ~ e 2 /h. Whether this value is equal to the 
above universal conductivity 4e 2 /wh, Eq. I|79(l . or the nu- 
merical coefficient is different, requires further study. Re- 
markably, the same value Eq. I|79|) has been founcliA*22i2£ 
for the conductance of a clean graphene sample in the 
considered geometry L y ^> L x . 

3. Cylindric geometry 

Let us take a Co-symmetric strip of a large transverse 
size L y and infinite in the x direction, and roll into a 
cylinder, preserving the chiral structure. Let us further 
assume a small but non-zero level width 5, as in Sec. llV Al 
If 5 is much larger than the level spacing in the square 
L y x L y , the system is effectively two-dimensional, and 
the consideration of Sec. IIV Al applies. Let us consider 
the opposite limit. One can then ask whether our result 
concerning the universal conductivity will be applicable 
in this quasi-lD geometry. Analyzing the derivation in 
Sec. llV Al it is not difficult to see that it breaks down: the 
momentum q y is now quantized and its shift therefore not 
allowed. Let us consider, however, an Aharonov-Bohm 
flux <!> piercing the cylinder, which amounts to introduc- 
ing the extra phase e J */*° in the periodic boundary con- 
ditions ($o is the flux quantum). Averaging over we 
restore the applicability of the consideration of Sec. llV Al 
so that 

4e 2 

(tr), = — . (95) 

7Tft 

Therefore, depending on the value of the Aharonov- 
Bohm flux, the conductivity can be either larger or 
smaller than this universal value, which is restored af- 
ter the averaging. A similar strong dependence of con- 
ductance of a clean graphene strip on the boundary 
conditions was found in Refsiii^. Our observation of 
the Aharonov-Bohm flux dependence of the conductivity 
seems also to be related to the known results on trans- 
port properties of disordered wires with chiral symme- 
try, namely their dependence on the parity of the num- 
ber of channels and the staggering in the hopping matrix 
elements 



V. CONCLUSIONS 

To summarize, we have studied electron transport 
properties of a disordered graphene layer. We have shown 
that the nature of disorder is of crucial importance for 
the behavior of the conductivity. Away from the half 
filling, the concentration dependence of conductivity is 



linear (with logarithmic corrections) for strong scatter- 
ers (unitary limit), while it is only logarithmic in the 
case of weak scatterers (Gaussian disorder). We have 
constructed a "phase diagram" showing which of these 
types of behavior should be expected for given micro- 
scopic parameters of the disorder. We have shown that 
the physically important case of charged impurities corre- 
sponds to the Gaussian-unitary "phase boundary" . The 
linear behavior of the conductivity that we have found for 
the case of strong scatterers agrees with the experimen- 
tal findings, 2,3 demonstrating that this kind of disorder 
is dominant in experimentally studied structures. 

At half filling, the conductivity is of the order of e 2 jh 
if the randomness preserves one of the chiral symmetries 
of the clean Hamiltonian; otherwise, the conductivity is 
strongly affected by localization effects. For the case of 
chiral disorder, the exact value of the conductivity still 
depends on the nature of the infrared cut-off, which may 
depend on the physical setup. We have analyzed in de- 
tail the situation when this cut-off is provided by the 
level width 8 or by the frequency w; in the first case the 
conductivity takes a universal value 4e 2 /irh, while in the 
second case it shows a more complex behavior. Whether 
the chiral disorder may indeed dominate in experimen- 
tally relevant structures, explaining the observed value 
of conductivity ~ e 2 /h remains an open question. From 
the theoretical point of view, further research directions 
extending our results include, in particular, the meso- 
scopic transport in a phase-coherent disordered sample 
and effects of interaction. 
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APPENDIX A: ONE-LOOP RG EQUATIONS 

A complete set of one-loop perturbative RG equations can be obtained by considering the diagrams of Fig. [5] with 
all possible disorder structures from Tabled An impurity line in those diagrams represents a sum of all possible types 
of disorder with the proper amplitude and corresponding matrices at the vertices. The RG equations for 9 disorder 
amplitudes [diagrams (b) , (c) , and (d) in Fig. |S] have the form 



2a (a + Po + 70 + a± + P±_ + 7_l + a z + Z + j z ) + 2a±a z + p±p z + 2~/ ± j z , (Ala) 

2{2a a z + o p z + 2 7o7z ), (Alb) 

-2a z (a +/3 +j -a±-/3±-i±+a z + l3 z + j x ) + 2a a± + (3 (3 ± + 2j a j±, (Ale) 

2[/3 (ao - 7o + + ol z - j z ) + a±(3 z + a z f3± + fixlo], (Aid) 

4(a /3 z + a z f3 + [3 ~/ + Pxlx + &7*)i ( Ale ) 

2[-P z (a ~ 7o - + a> z - j z ) + a {3 ± + a±(3 + 0xlz], (Alf) 

27o(«o - A) +7o + "J. - Px + Ix + ol z - /3 Z + j z ) + 2a ± j z + 2a z -f ± + O P±, (Alg) 

4a o7z + 4a z7o + (3 2 + + 01 , (Alh) 

-27^(«o - ax + a z ~ Po + Px - Pz + 7o - 7_l + 7a) + 2a o7 _L + 2a ±7o + PxPz- (All) 



dlogi 

The RG equation for the energy [diagram (a) in Fig. [S] reads 
ds 



dlogL 



e(l + a + Po + 7q + a± + p± + -y± + a z + p z + lz ). (Alj) 



For brevity, in this Appendix we omit tildes which distinguish running parameters from their initial values in the 
main text. 

In various particular cases, when only some subset of disorder structures is present, these equations reduce to the 
corresponding form known in the literature. The cases of Co-chiral (a±, fix, ~f±) and C z -chiral (a±, Po, p z , 71) 
disorder are considered in Sec. II V Al If the disorder is proportional to the T3 matrix (ax, 7oj lz), the Hamiltonian 
decouples in two 2x2 blocks, which have the structure of the model with random mass (7 Z ), scalar (70), and vector 
(o>±) potential analyzed in Ref. l60t The RG equations for the random mass problem were also given in Refs. I59ll62l 
for the random potential in Ref. l6ll 

If the system possesses a time- reversal invariance (To), only the couplings ato, 0±, P z , 7j_, and 7 Z survive, which 
is the case considered in Ref. |3(1 Taking into account the difference between our RG scheme and that of Ref. |3(] 
(where the velocity is renormalized whereas the energy is not), we have checked that RG equations of Ref. |3(] are 
reproduced from the complete set l|Alf) if a number of assumptions concerning the hierarchy of the disorder couplings 
(q > V a o|2/3 z -Px\, v / "o|27 z -7j_| > p z ,Px,7z,7x > \20 x - 0x\, \^ z - Jx\) are made. 

APPENDIX B: IMPURITY-INDUCED CORRECTIONS TO THE DOS IN THE UNITARY LIMIT 

In this Appendix we calculate the density of states in the presence of infinitely strong impurities (unitary limit) 
up to the second order in their concentration n; mp . The contribution of the first order in ni mp is determined by the 
diagram (Fig. |Bk) containing a single T-matrbi§*i 

(1) 4 Tm 77A 2 f d 2 k e*+vtk 2 _ rjA 2 



SP {£) 7r Im 2[elog(A/|e|) + z^| £ |/2]7 (2tt)2 [{e + i0) 2 _ v 2 k2] * 2 ™ 2 |e| log 2 (A/| £ |) ' (B1) 

This result corresponds to Eq. JSSJ. Obviously, the first-order contribution to the density of states is correctly taken 
into account by the self-consistent unitary approximation. 
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FIG. 8: Diagrams for (a) the first-order correction to the density of states and (b) the second-order correction to the partition 
function. 



The problem becomes more complicated when one looks for the second-order contribution. The calculations are 
greatly simplified in the coordinate representation and for Matsubara energies. The Green function and the T-matrix 
have the form 



G (ie,r) 



T(ie) 



2irv 2 



2nv 2 



i6elog(A/e)' 



1, long-range; 
4, short-range. 



(B2) 
(B3) 



The T-matrix has different values in the limits of long- and short-range potential disorder. Eq. i|B2|l for the Green 
function applies for not too short distance. One has to cut the real-space integrals at r ~ Vq/ A. 

We are going to express the density of states in terms of the partition function per unit area. The contribution to 
this quantity of the second order in rj is given by the diagrams Fig. [Bp 



Z 2 {ie) = nf mp Tr ]T — / d 2 r[G(r)G(-r)r = -2nf mp / d 2 r\og 1 
m=l m { 

The correction to the density of states can be represented in the form 
Sp^(e) 



1 Im dZ 2 

7r d(ie) 



ie— >e-\-iO 



Te 
2nv 2 



dz z K 2 (z) 



K ? - - Kl - 



b 2 \og'(A/e) + K 2 (z)-K 2 (z) 



(B4) 
(B5) 



ie— *e+zO 



For C — 61og(A/e) 3> 1 we split the integral over z into two parts and observe that it is dominated by the domain 
z>l/C: 



dz z K 2 (z) 



o £ 2 + K\iz) - K$(z) 



l/C 



dz z 



C 2 



dz z K\{z) 



l/C 



C 2 ' 



(B6) 



Substituting this in Eq. I|B5|) and performing the analytical continuation, we finally arrive at 



5p 2 (s) = 8v n- hn 



i loglog(A/e) 



_e 3 b 2 log 2 (A/e) 
where au(s) is determined by Eq. I|60|l . 



ie— *e+iO 



Wj4p loglog(A/|e|) 
6 2 |e| 3 log 3 (A/|e|) 



A 



= -2p(e)afj(e) log -- log log — , (B7) 
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